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Abstract. The investigation of the statistical properties of maps of line centroids has been used for almost 50 years, but there 
is still no general agreement on their interpretation. 

Aims: We try to quantify which properties of underlying turbulent velocity fields can be derived from centroid velocity maps, 
and we test conditions under which the scaling behaviour of the centroid velocities matches the scaling of the three-dimensional 
velocity field. 

Methods: Using fractal cloud models we study systematically the relation between three-dimensional density and velocity fields 
and the statistical properties of the produced line centroid maps. We put special attention to cases with large density fluctuations 
resembling supersonic interstellar turbulence. Starting from the A- variance analysis we derive a new tool to compute the scaling 
behaviour of the three-dimensional velocity field from observed intensity and centroid velocity maps. 

Results: We provide two criteria to decide whether the information from the centroid velocities directly reflects the properties 
of the underlying velocity field. Applying these criteria allows to understand the different results found so far in the literature 
on the interpretation of the statistics of velocity centroids. The new iteration scheme can be used to derive the three-dimensional 
velocity scaling from centroid velocity maps for arbitrary density and velocity fields, but it requires an accurate knowledge of 
the average density of the considered interstellar cloud. 
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^ 1 . Introduction direct way to do so. Observations of the profiles of atomic or 

molecular lines from interstellar clouds allow to deduce infor- 

. ^ Understanding the role and nature of interstellar turbulence has mation on the line-of-sight velocity structure of the clouds. The 

^ , been the subject of intensive studies for half a century now problem of recovering of the velocity information from lines is 

^ ' but st ill remains open in many aspects ("cf. lElmegreen & Scalol far from being straightforward. Even in the most simple case 

l2004l) . Major questions concern the mechanisms by which tur- of thermally excited optically thin lines from an isothermal 

bulent motions are driven and the role of the strong compress- medium the line profiles originate from a convolution of the 

ibility of the interstellar medium on the structure of the tur- density structure p depending on the sky coordinates x = (a, 6) 

bulent energy cascade. Both aspects are directly reflected on and the line-of-sight coordinate z with the velocity structure 

the spectrum of velocity fluctuations in the turbulent motion. v z (x,z): 
It is frequently claimed, that driving mechanisms should create 

dominant motions at the corresponding scales, and the power v ) oc I dzp(x z)<f>(v - v z (x z)) (1) 

spectrum of velocities in the turbulent cascade is known to J 

change from a PQk\) oc |*|-»/3 Kolmogorov spectrum for an in- ^ ^ of nam)w ^ ^ ^ profile y ( ^ z)) cm bg 

compressible medium to a P(\k\) oc \lc\-* spectrum of Burger's approximated by a &iunctioQ . The re are seve ral complemen- 

turbulen ce in a highly co rrirjressible medium dominated by tary wayf , tQ ufje Ms information ( c f jLazariarll20ol) . Here, we 

shocks < Chappell& Scalol U223- However, numerical simu- restrict ourselves to cen troids, the first moment of the lines, but 

lations show often a different behaviour (see | Cho & Lazarian[ ^ centmids m nQ direct map of ^ vdodty stmc _ 

2005) which makes the issue of the observed spectrum very t 

intriguing. Models for the density structure and the relation between 

To support the theoretical understanding of the interstellar density and velocity structure are needed to deduce the latter 

turbulence it is thus essential to actually measure the velocity from the profiles I(x, v). This is straight-forward for simple ge- 

structure in the interstellar medium. Unfortunately, there is no ometries like spherical clouds or thin disks but extremely diffi- 
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cult for filamentary turbulent cloud structures showing varying 
substructures on all spatial scales. Hydrodynamic or magneto- 
hydrodynamic numerical simulations can be used as physically 
justified models for turbulent interstellar clouds within a lim- 
ited dynamic range. The nature of these can be described, how- 
ever, only in terms of statistical measures. Fractal cloud mod- 
els provide a reasonable phenomenological description of the 
clouds. We focus on measures for the spatial scaling of the 
velocity structure. The ultimate goal is to derive the three- 
dimensional (3-D) power spectrum of velocity fluctuations. 

A recovery of 3-D information from the available 2-D 
data requires, in general, an inversion, which may result in 
substantial noise in the inverted data. For deriving the tur- 
bulence statistics we can, however, use its symmetries. Here, 
we restrict ourselves to statistically isotropic turbulence. The 
derivation of properties of anisotropic, but axisymme t ric tur - 
bulence from observations was discussed by iLazariaij (| 1 995 ) . 
Anisotropics can be due to the magnetic fields faig don , 
ll984llZank & MatfhaeusL Il992t Idoldreich & Shridhaiill995l) . 
How ever, if, as both theo r y and numerics sugg ests (see 
Idoldreich & Shridharl Il995t ICho & Lazarianl l2003l) . the en- 
ergy spectrum is dominated by fluctuations perpendicular to the 
local direction of magnetic field, the effe cts of anisotropy on 
the observed spectra can be neglected (Esa uivel etaHl2003l) . 

In order to derive the isotropic power spectrum we will use 
an auxiliary quantity, the A-variance spectrum, because of prac- 
tical advantages when measuring the velocity scaling in ob- 
served data. Moreover, we restrict the analysis here to the first 
moments of the lines, the centroid velocity, as the most obvious 
tracer to measure the velocity structure in an interstellar cloud. 

Maps of observed line centroids have been systematically 
studied to obtain the scaling behaviour of centroid velocity dif- 



Munch, 


1958; 


iKlginer^Dickmgril 


Il985 


1994; Lis etal., 


1996; Miesch et al., 


1999) 



Mies ch & Balh 
However, there is 
still no agreement on the theoretical relation between the ob- 
served scaling behaviour of the centroid velocities and the scal- 
ing behaviour of the underlying turbulent velocity structure. 
Although it was clear from the very beginning that density 
structure can influence the line centroids, up to the recent past 
there was no criterion to estimate to quantitative effect of den- 
sity. 

Investigating hydrodyna mic turbulence simulations 
Ossenkop f & Mac Low! J2002I) found that the centroid maps 
show approximately the same Hurst index, i.e. the same 
relative variation across a given scale, as the underlying 3-D 
velocity structure. This means that the power spectral index in 
their centroid maps was reduced by one compared to the power 
spectral i ndex in the 3-D velocity structur e. Studies of fractal 
clouds by Miville-Deschenes et al. (2003a) showed in contrast 
that their centroid maps show th e same power spectra l index 
as the 3-D velocity structure 1 lLazarian & Esa uivel (2003) 
provided an analytical treatment of the centroid statistics 



introducing a new more robust definition of velocity centroids, 
formulated a criterion when the centroids represents the veloc- 
ity statistics, but this publication did not cover the parameter 
space to be fairly compared with previous studies. T he prob- 
lem w as further elaborated in a subsequent study by iLevrierl 
( 2004), who pointed out that the statistical t reatment presented 
in the form of structure functions by lLazarian & EsauiveH 
(2003) may have some advantages if rewritten in terms of 
correlation functions. Assuming that the fluctuations are small 
compared to the mean density he obtained analytic expressions 
for correlation functions of centroids. Combi n ing str ucture and 
correlation functions Esau ivel & Lazariai] J2005I) provided 
a detailed study of centroid velocities for data obtained 
through compressible MHD simulations. Here, we compare 
the different centroid definitions and test their outcome for a 
set of fractal cloud models. 

Using the A-variance analysis of the centroid maps, we 
show that it is in principle applicable to derive the velocity 
power spectrum from observed centroid maps but that the reli- 
ability of this derivation depends critically on individual turbu- 
lence parameters. The centroid maps reflect the actual velocity 
distribution only in a medium with an average density which is 
large compared to the density dispersion. Here, the A-variance 
analysis provides a direct measure for the power spectral in- 
dex of the velocity structure. Only when applied in an iterative 
process with an a priori knowledge of the average density, the 
analysis of centroid maps allows to approximate the velocity 
structure in the general case. The approximation is better as 
steeper the velocity spectrum is and as better the average den- 
sity is known. 

In Sect. 2 we shortly repeat the formalism used to describe 
the velocity centroids, discuss the properties of the test data 
sets, and the ways to measure their spatial scaling behaviour 
in terms of the A-variance. In Sect. 3 we perform the analysis 
of the centroid maps using the A-variance, compare the results 
with the original test data and derive criteria when the centroid 
maps can be used to measure directly the three-dimensional 
velocity structure. In Sect. 4 we propose an iterative method 
to derive the power spectrum of the velocity structure from the 
centroid maps in cases without a direct matching. Sect. 5 gives 
a summary with respect to the interpretation of observed data. 

2. The starting point 

2. 1. Definition of centroid velocities 

For the fluctuating density and velocity fields in a cloud we can 
always write 



p(x, z) = po + <5p 
v(x, z) = vq + Sv 



(2) 



where po and vo are averages over the whole cloud, and dp and 
5v denote the variations across the cloud 2 . 



1 When dealing with projected quantities one has to carefully dis- 
tinguish correlation functions and power spectra. When a power-law 
approximation is good for both of them, the spectral index of corre- 
lation functions gets steeper by one due to projection, while the 2-D 



projected power spectrum retains the spectral index of the underlying 
3-D spectrum. 

2 From here on we drop the index z in the notation for the line-of- 
sight component of the velocity because we consider only this com- 
ponent. 
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When we assume that the emissivity is proportional to the 
density of the cloud, the line intensity I(x,v) at velocity v is a 
measure for the total column density of emitters with this ve- 
locity at a given line-of-sight x. This condition is violated for 
optically thick lines or media with strongly varying tempera- 
tures but it is e.g. well fulfilled for the [CII] emission from the 
cold neutral medium or the HI emission from the warm neu- 
tral medium. The effect of self-absorption will be quantified in 
a subsequent paper. For constant emissivity the integrated line 
intensity is 

/intOO = J dv I(x,v) 
= X J dzp(x,z) 



(3) 



where X is the proportionality factor from Eq. translating 
the column density into a line intensity. 

There are two different centroid definitions in common use. 
Ordinary centroid velocities, also known as normalised cen- 
troids, are obtained as 

(x) = J dvvl(x,v) J dvl(x,v) (4) 
= J dz (v'o + Sv)(p + 6p)j j dz(p + dp) 

Unfortunately, this definition implies a complex combina- 
tion of density and velocity fluctuations which makes it im- 
possible to disentangle the influence from both structures in 
the general case. Only in case of very s mall fluctuation s, a lin- 
earisation technique can be developed dLevrierl.|2004l) . A bet- 
ter separation of density and velocity fluctuations in the cen- 
troids is obtained when we apply the definit ion of weighted, i.e. 
unno rmalised, centroids as proposed by lLazarian & Esauivel 
(200§ 3 



v c (x) = l/xj dvvl(x,v) 
= VoPoZtot + Po J dz Sv 



jdz6 P + j 



dz 5p5v 



(5) 



where z tot is the total thickness of the cloud. In this definition 
the centroids do not have the dimension of a velocity but of 
velocity times column density. For a better comparison with 
the ordinary centroid velocities it is useful to normalise the 
weighted centroids by the average column density poZtot> but we 
omit this factor in the following to keep the equations shorter. 
The constant factor would not change any of our conclusions 
on the scaling behaviour of the velocity structure. 

We see that even in this definition the centroid velocities are 
not simply determined by the projected velocities vq + f dz dv 

3 In contrast to the original definition we have not included the con- 
stant factor X in the centroid definition so that the weighted centroids 
have the dimension of a velocity times column density here instead 
of velocity times intensity. This keeps the equations in the following 
sections somewhat shorter. 



but also by two terms reflecting the density variations. The con- 
tribution from the projected density variations J dz dp can be 
easily obtained from the integrated line profiles and it can be 
eliminated by selecting a velocity scale with vrj = 0. However, 
the term containing the product of the fluctuations in the den- 
sity and the velocity structure cannot be measured separately. 

The scaling behaviour of the centroid velocities depends 
on the combination of density and velocity variations along the 
line of sight, which cannot be retrieved directly. The relative 
contribution of the simple projection of the velocity structure 
and the density variations 5p across the line of sight depends 
on the ratio between the density fluctuations 5p and the average 
density po. Eq. (|5jl thus shows already that the ratio between the 
density dispersion cr p and the average density po is a critical 
parameter for the relation between the 3-D velocity scaling and 
the centroid scaling. 

2.2. Test data sets 

To study the general ability of different methods to extract the 
underlying velocity structure from observed centroid velocities 
we construct well defined test data sets for the density and ve- 
locity structure which are used to study the translation of their 
scaling properties into centroid properties. 

Interstellar cloud observations often reveal self- similar 
scaling properties (e.g. Falgarone et all 119 95: Combes, 2000) 
corresponding to power-law power spectra of the intensity 
distribution. Such intensity maps can be approximately mod- 
elled by fractional Brownian motion (fB m) structures (see e.g. 
IStutzki et all Il998t iBensch et all 1200 lb . They are defined by 
the single number /? determining the exponent of the power 
spectrum, P(\k\) oc \k\~P. The phases of the Fourier spectrum 
are random. 

Thus fBm's represent one of the simplest possible repre- 
sentations of interstellar cloud structures still allowing a pa- 
rameter study in terms of the spectral index /3 which deter- 
mines the actual appearance of the structures. fBm's can be 
defined in arbitrary dimensions and we use their essential prop- 
erty that the projection of an fBm to lower d imensions re s ults in 
a new fBm with the sam e spectral index llSmtzkiet .all 1 1998: 
iBrunt & Mac Low! I2004I) . 4 Thus the spectral index measured 
for the column density directly reflects the index of the three- 
dimensional density structure. 

Measured spectral indices for the column den- 
sity structure of interst el lar cl ou ds range from 2.0 t o 
3.7 felmegreen & ScaloL l2004 iFalearone et ail l2004 . 
Observations of large molecular clouds and molecular cloud 
complexes and HI absorption line studies provided ty pical 
value s between 2.4 and 2 . 9 (e.g .lstenholmil l984[lLan ger et all 
1993[ iDeshpande eTall l200ot IBenschetall l200lt iHubeil 
2002; Padoa netall l2003h . whereas IBenschetall foOOlb 



found indications for somewhat larger indices at the scales of 
cloud cores. Observations of the w arm atomic gas prov ided 
typical values between 3.3 (e.g. Stani mirovic & Lazariarl 



4 As discussed by Stutzki et al. 1 1998) it is easy to show that this 
fBm property violates the often used hypothesis th at the fract al di- 
mension decreases by one in projection (Peitgen & SaupeL ll988). 
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2001) and 3.6 dMiville-Deschenes et all l2003bl) with some 
indications for an even broader range from 8/3 to 11/3 in the 
LMC Elmegreen2001. Due to a lack of direct measurements, 
as discussed in the introduction, the index range of the 
velocity structure is sti ll hardly known. MHD simulations by 
ICho & Lazarianl J2003) indicate that i t should be close to the 
Kolmo goro v value of 1 1/3. I n cont rast lOssenkopf & Mac Low! 
J2002I) and brunt & Heverl J2002h obtained velocity spectral 
indices close to four from observations of the Polaris Flare 
molecular cloud and of molecular clouds in the FCRAO 
survey of the Outer Galaxy, respectively, consistent with the 
properties of a shock-dominated medium. In these cases, the 
velocity spectrum was always steeper than column density 
spectrum. Here, we do not aim at reproducing the exact 
combination of spectral indices for any particular interstellar 
cloud, but want to study the general behaviour covering the 
ful l range of spectral ind i ces o bserved so far. 

Esau ivel & Lazarianl J20 05) demonstrated that the centroid 
structure function shows a qualitatively different behaviour for 
spectra with an index above and below 3.0 (steep and shallow 
spectra). Unfortunately, the observational data do not rule out 
either of the two types. Thus we focus on two test data sets: 
fBm's with a spectral index of 3.7 representing steep spectra 
and with an index of 2.6 representing a shallow behaviour, re- 
spectively. They sample both regimes and are close to some 
observed values for the velocity and density structure. We have 
studied a much larger parameter range covering spectral indices 
between 2.0 and 4.0 but with the four possible mutual combi- 
nations of the two mentioned spectral indices all major effects 
are covered so that we restrict ourselves to these cases for all 
examples given in the following. 

In Fig. [T] we give a visual impression for the difference in 
the actual projected structure between fBm's of different spec- 
tral index. The spectral index basically determines the relative 
contribution of structures on different size scales. The ffim with 
an index of 2.6 shows a large amount of small scale clumps and 
filaments, whereas the ffim with [5 — 3.7 consists basically of 
one peak with fragmented boundaries. 

The figure also reveals a general problem of ffim's when 
interpreted as density structure. They show negative val- 
ues, ffim's have on the average a Gaussian probability dis- 
tribution with vanishing mean so that negative values can 
only be avoided when adding a large constant density offset. 
However, in this way we drastically change the ratio <r p /po 
for the data set. Another method to create a density distri- 
bution containing only positive values i s to square or expo - 
nentiate the original ffim as pro posed bv lStutzki et alJ dl998). 
Miville-Deschenes et alJd2003al) have claimed that exponentia- 
tion, p exp = po exp(apfB m ), does not affect the power-spectrum, 
but it is mathematically obvious, that it can potentially destroy 
the power-law scaling. Thus we have tested the impact of expo- 
nentiation for different spectral indices /? and different factors a 
translating the standard deviation of the ffim into the logarith- 
mic standard deviation of the new density structure. The result 
is shown in Fig.Elfor an ffim w ith (3 — 4 as used in Fig. 1 1 of 
iMiville-D eschene s et alJ d2003al) and for an ffim with /3 = 2.6 
in terms of A-variance spectra (see Sect. 12. 3> . It is obvious that 
for narrow distributions, the distortion of the original spectrum 




20 40 60 80 100 



Fig. 1. Projected maps of ffim structures with spectral indices 
/3 = 2.6 (upper plot) and f3 - 3.7 (lower plot). Both data sets 
use the same random phases leading to the apparent similarity 
of the overall distribution in this example. 

by exponentiation is small, as the exponentiation is then close 
to a linear transformation. In general, we have to acknowledge, 
however, considerable distortions of the spectrum by the expo- 
nentiation. When creating a very wide density distribution from 
the /3 = 2.6 ffim we even find a completely different scaling 
behaviour resembling rather a structure wi th /3 — 0. The ex- 
ample from lMiville-Deschenes e t al. (2003a) corresponds ap- 
proximately to the (3 — 4, <r\ og = 0.3 case shown in Fig. |2] 
Here, the deviation from the original spectrum is small so that 
it was not detectable. Moreover, we have found that the A- 
variance reacts much more sensitive to the exponentiation than 
the azimuthally averaged power spectra. Only for very wide 
distributions and low spectral indices the azimuthally averaged 
power spectra show similar noticeable deviations. In general 
we have to conclude that exponentiation results in a change 
of the scaling properties. Consequently, non-linear transforma- 
tions are not well suited to produce well defined test data for 
the density structure. We will stick to the simple approach of 
adding a constant to the ffim and ignoring the remaining neg- 
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2.3. The A-variance 
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Fig. 2. A-variance spectra of the projected structure of expo- 
nentiated ffim's compared to the spectrum of the original fBm. 
The upper plot represents (3 — 4, the lower plot ft = 2.6. The 
different lines indicate different stretching factors a resulting 
in different logarithmic widths of the distributions. The aver- 
age logarithmic density is taken to be 2.0 in all cases. 



ative values for the construction of the density structure. The 
implications of this approach are quantified in detail in Sect. 
13.11 In contrast to the density structure which has to be positive 
defined and necessarily a non zero mean po, the velocity struc- 
ture can directly use fBm's thus guaranteeing a zero value for 
vo so that the simplifications discussed above apply. 

When using independent fBm's to represent both the den- 
sity and the velocity structure of interstellar clouds we ne- 
glect, however, the interrelation of both quantities in the 
interstellar medium determined by the hydrodynamic equa- 
tions, especially by the Poisson equation. Comparisons with 
magneto-hydrodynamic simulations by lEsauivel & Laz arian 
(2005) have shown, however, that the cross-correlation be- 
tween the density and velocity fields has a negligible effect on 
the centroid velocities so that we can neglect its impact here. 
We will further discuss the influence of cross-correlations be- 
tween density and velocity structure on different observational 
parameters in a subsequent paper. 



The A-variance analysis was introd uced bvlStutzki et alJ|l 998) 
and improved and ex tended by iBensch et alJ (1200 lh and 
Ossenkopf et al. (2005). Here, we repeat only those definitions 
which are essential for the centroid analysis. 

The A-variance in a structure fix) is computed by filtering 
the data set with a spherically symmetric normalised wavelet 
of characteristic size /, consisting of a positive inner part and 
a negative annulus, and computi ng the variance of the fil- 
tered map. lOssenkopf et alJ ( l2005l) have tested various wavelet 
shapes, but their mutual differences are not significant for the 
analysis performed here so t hat we stick to the ordinary French 
hat filter from Stut zki et alJ ft 998). The A-variance is then the 
variance of the filtered map, as a function of the filter size, 
given by 



aim = ((fix) * ©,(*)) : 



(6) 



where the symbol * stands for a convolution, Q, describes the 
filter wavelet and the average is taken over the whole data set. If 
I is the average distance between two points in the core and the 
annulus in the filter, the A-variance spectrum <x?(/) measures 
the amount of structure on the given scale /. 

The A-variance is related to the power spectrum of a struc- 
ture P(k) by 



ol(/) = j P{k) 0/i*d 



d"k 



(7) 



where Q, is the Fourier transform of the filter function with 
the size I and k denotes the spatial frequency or wavenumber. 
In case of isotropic structures the power spectrum is spherically 
symmetric, P(k) = P(\k\). This is also the case for the Fourier 
transformed filter function as long as it is spherically symmet- 
ric in the spatial domain. The power spectrum is given by the 
Fourier transform of the autocorrelation function 

A(l) = (f(x) f(x + t)) r (8) 

For power-law power spectra Istutzki et al. I il998t> showed 
that for 2-D structures in the interval of spectral indices /3 be- 
tween and 6, the A-variance spectrum is as well a power law 
with the exponent a = /3 - 2. In three dimensions the range 
is extended to < f3 < 7 and the exponent is a = y8 — 3. 
Equivalent slopes are obtained locally in case of non-power- 
law power spectra. However, in this case there is no analytic 
relation for the normalisation factor of the A-variance spectrum 
so that it can only be obtained by numeric integration. 

Thus the A-variance is basically a very robust method to 
evaluate the power spectrum of a structure. The advantages of 
the A-variance compared to the direct computation of the power 
spectrum result from the smooth filter shape which provides a 
very robust way for an angular average independent from grid- 
ding effec ts, and from the ins ensitivity to edge effects as dis- 
cussed bv lBensch et alJ ( l200ll) . A possible disadvantage is the 
implicit radial averaging, which does not allow to search for 
signatures of anisotropy still contained in the two-dimensional 
power spectrum P(k ) . Such an anisotropy was considered by 
Esquivel & Laza rian ( 2005J) but is irrelevant for our studies. 
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2.4. Comparing A-variance and structure function 

lLazarian & Esauivel d2003l) and lEsauivel & Lazarianl d2005h 
used the (second order) structure function instead of the A- 
variance to characterise the scaling of velocity centroids. The 
structure function is related as well to the autocorrelation func - 
tion, D(l) = 2[A(0)-A(/)], (see e.g. iMiesch & Ballvi fl994l) . 
With the power spectrum being the Fourier transform of the 
autocorrelation function we also have a trivial relation between 
structure functions and power spectra. 

For stru c tures with a power-law power spectrum 
IStutzki et al.l Jl998h studied analytically the relation be- 
tween the power spectrum, the autocorrelation function and 
the A-variance. They find in the range of spectral indices 
3 < p < 5 in 3-D and for 2 < /3 < 4 in 2-D, and in the limit 
of infinitely large data sets, power-law stru c ture f unctions. 
Using the notation of lLazarian & Pogosvanl (l2000l) . this is 
the range of steep spectra. Here, the spectral index of the 
structure function agrees with the index of the A-variance 
spectra discussed above. In the range of shallow spectra 
with lower power spectral indices, 0<y8<3or0<y8<2 
respectively, the autocorrelation function is a power law 5 so 
that the structure function must deviate from a power law 
behaviour. The structure function is always increasing with lag 
towards the maximum given by twice the total variance of the 
structure cri = A(0). 

For MHD simulations producing basically steep velocity 
spectra but with signifi c ant deviations from pure power laws 
lOssenkopf & Mac Low! d2002h compared the centroid veloc- 
ity structure function with the A-variance of the centroid map 
and showed that both give a similar scaling behaviour, having 
comparable slopes within a large part of the spectrum., The A- 
variance, however, is advantageous with respect to the detec- 
tion of pronounced scales in the map and is more robust with 
respect to observational artifacts. Altogether, the A-variance 
seems to be somewhat better suited to determine the exponent 
of the power spectrum, as it shows itself a larger range of a 
power-law behaviour and it is more stable with respect to ob- 
servational restrictions. 

On the other hand, lEsauivel & Lazariai] J2005I) demon- 
strated that the structure function of centroid velocities can be 
analytically understood with respect to its composition from 
density and velocity fluctuations. This represents a clear ad- 
vantage relative to the A-variance. Thus we have actually per- 
formed all tests of the centroid structures reported here both 
with the A-variance analysis and with the structure function. As 
a surprising result we find very little differences in the general 
behaviour. Therefore, we concentrate in the following analysis 
on the A-variance spectra discussing the differences in compar- 
ison to the structure functions only in Sect. 13. 51 

2.5. Projection effects 

The relation between a 3-D structure and projected 2-D 
maps, obtained by the integration along the line of sight, has 

5 For spatial separations corresponding to wavenumbers smaller 
than the cut-off wavenumber given by the finite sampling of any sys- 
tem. 



bee n studied in detail both in terms of the A-va riance (e.g. 
IStutzki et all 1 19981 iMac Low & Ossenkopfl l2000l) and of the 
structure function ( Esau ivel & Laz arian. 2005). Aprojection of 
the density structure p(x,z) is inherently performed when ob- 
serving the intensity map Ii nt (x) of an optically thin tracer in a 
medium of constant excitation temperature (Eq.^. 

The projection effect on the A-variance spectrum can be 
easily understood by realizing that the A-variance is basically a 
robust method to deduce the power spectrum. In Fourier space, 
projection corresponds to the selection of the zero-frequency 
component in the considered direction. For isotropic structures 
the power spectral indices of projected maps in any direction 
agree with the spectral index of the 3-D structure. This is ful- 
filled by definition for the fflm structures used here for testing. 
Thus the local slope of the power spectrum j3 is retained, and 
all components with non-zero spatial frequencies in the consid- 
ered direction are dropped. Because the A-variance is obtained 
by convolving this power spectrum with the Fourier transform 
of either a 3-D or a 2-D wavelet, the resulting spectrum has a 
local slope o^d = /3 - 3 or 020 =0 — 2, respectively. The mu- 
tual translation is straight forward. The exponent of the power 
spectrum is retained on projection, while the index of the A- 
variance 6 is increased by one. This has been confirmed in the 
application of the A-variance analysis to the 3-D density struc- 
ture of (magneto- )hydrodynamic simulations and their projec- 
tion onto maps by M ac Low & Ossenkopil 1E0OO). 

For power-law power spectra, the translation of the ampli- 
tudes can also be performed analy tically follow i ng the formal- 
ism provided in the Appendix of IStutzki et alJ dl998t) . As an 
approximation we can also use the simple empirical relation 

< 3D (0 = < 2D (0 x jL. x 1.97 exp(-^) (9) 

which is accurate within a few percent for power spectral in- 
dices (3 between 1 and 4 and cube sizes of at least 32 3 pixels. 
Even for sufficiently smooth, but non-power-law A-variance 
spectra Eq. l|9} can be applied by using an index f3(l) derived 
from the local slope. 

A general problem is, however, the actual loss of infor- 
mation by projection. There is no way to recover the Fourier 
amplitudes which are dropped by the projection. Thus the re- 
translation from the 2-D A-variance spectrum into the corre- 
sponding 3-D spec trum i s only possib le by assuming isotropy. 
IMac Low & Ossenkopl d2000l) and lOssenkopf & Mac Lowl 
(2002) studied the degree of anisotropy in hydrodynamic and 
magneto-hydrodynamic simulations by comparing 2-D and 3- 
D A-variance spectra and found that the assumption is clearly 
violated for simulations with strong magnetic fields but reason- 
ably justified for most other simulations. 

Fig- demonstrates the influence of the projection effects 
on the A-variance spectra of two fBm's. The upper graph rep- 
resents an fBm structure with a shallow index yS = 2.6 and the 
lower graph a steep spectrum with f3 = 3.7. The A-variance 
spectra measured in 3-D and for the projected structure follow 
almost exactly the theoretical power-law relation with the ex- 
ponents a = yS — 3 or a = ft— 2, respectively. The triangles stand 

6 The same applies to the structure function, but in a limited spectral 
range. 



Ossenkopf et al.: The statistics of centroid velocities 



7 




rescaled 3-D A— variance A 
3-D A-variance x 
A-variance of projection 



lag [pixel] 




lag [pixel] 

Fig. 3. A-variance spectra determined in 3-D structures and 
their projection together with the translation of the 3-D A- 
variance spectrum into the corresponding 2-D spectrum using 
Eq. (|9}- For the upper plot an fBm with f3 = 2.6 shifted by 
2cr and truncated at zero was used. The lower plot shows the 
result for an fBm with p 1 = 3.7. The solid lines mark slopes 
corresponding to the spectral indices /3 - 3 and /3 - 2. 

for the results from the A-variance computed in 3-D and trans- 
lated into a 2-D spectrum using Eq. (|9j. We find an excellent 
agreement with the spectra obtained directly from the projected 
maps. 

Beyond the plotted range the A-variance spectra show a 
turn-over at about half of the total size of the simulated cube 
arising from the lack of larger structures due t o the periodic- 
ity co ndition in the construction of the data (see lBensch et all 
2001). Because of the lacking significance at large lags, the 
spectra are only computed up to lags of about a third of the 
cube size. 



3. Centroid composition effects 

Taking their relation to the autocorrelation function both the 
A-variance and the structure function of velocity centroid maps 
will be given by averages of the products v c (x)v c (x + 1) (see Eq. 
IS). Using the decomposition of the velocity centroids in Eq. (jSJi 
and assuming a zero average velocity vq we see that four terms 



characterise the scaling: 

+ | J" dz 6p(x, z)6v(x, z) x J dz 6p(x + I, z)5v(x + I, z)j 

+po | J dz 5v(x, z) x J dz 5p(x + /, z)6v(x + l,z)j 

+po ^ dz Sp(x,z)6v(x,z) x J dz 5v(x + l,z)j (10) 

The first term is the autocorrelation function of the projected 
velocity fluctuations. If this term dominates, the scaling be- 
haviour of the centroid velocities reflects exactly the scaling 
behaviour of the velocity structure. In this case it is easy to de- 
duce the properties of the velocity structure from an observed 
map of centroids. We find the simple projection of the velocity 
structure onto a 2-D map like in the case of the column density 
map reflecting the 3-D density structure. The second term de- 
scribes a combination of the fluctuations of the density and the 
velocity structure. The term also contains the mutual correla- 
tion between density and velocity fluctuations along the line of 
sight. The third and fourth terms quantify the cross-correlation 
between velocity fluctuations at one point and density fluctua- 
tions at another point. In case of isotropic media both terms are 
identical. They should statistically vanish in case of indepen- 
dent density and velocity structures, but some remainders due 
to accidental cross-correlations are expected for any particular 
realisation. A similar decomposi tion in terms o f the st ructure 
function was provided bv lLazarian & E sauivel (2003). From 
the decomposition in Eq. JlOb we see that the ratio between the 
average density and the density fluctuations should provide a 
criterion whether the centroid map is a good measure for the 
scaling of the velocity field. Thus we test in the following the 
composition of centroid velocity maps from fBm structures ad- 
justing their parameters in such a way that the full range of 
observed spectral indices for the density and velocity structure 
in the interstellar medium is covered. 

3. 1. The density zero level 

A major problem with the artificial simulation of density struc- 
tures is the mutual incompatibility of Gaussian fluctuations and 
strictly positive values for the density. As discussed in Sect. l2.2l 
fBm structures always show a Gaussian distribution of values 
and the analytic express i ons f or the velocity centroids derived 
by|Eazarian & Esquivell (2003) are also based on the assump- 
tion of Gaussian fluctuations. However, as long as the average 
of a Gaussian distribution is not large compared its dispersion, 
negative values are unavoidable for sufficiently large samples. 

A common way to create positive densities is to add a 
constant density until t he minimum value in the cube falls a t 
zero jMiville-Deschenes et all l2003at lEsauivel et ail l2003h . 
A major drawback of this method is, however, that the min- 
imum value of a Gaussian distribution depends on the ex- 
act realization of the random numbers used to generate the 
distribution and it is very sensitive to the size of the data 
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Fig. 4. A- variance spectra of 3-D density structures obtained by 
shift-and-truncate from an fBm with /3 = 2.6 and different trun- 
cation levels. The spectrum for the density cut at 3<x is partially 
indistinguishable from the original spectrum. 

cube. Thus the added value, then providing the average den- 
sity po, may significantly vary from simulation to simula- 
tion. By renormal ising the average densi ty to unity as pro- 
pose d by Miv ille-Deschenes et alJ J2003ah and lEsauivel et alJ 
( 2003p the variation is only transferred to the standard devia- 
tion of the density distribution because the ratio between stan- 
dard deviation and mean is retained. Mor eover, the approach 
results typically in <x p /(p) < 0.3 (Miville-Desch enes et all 
2003a). Such value s are in contr adiction to many observa- 
tional data (see e.g. Ijenkinsl 120041) . Density fluctuations with 
dplp ~ 1 are expected for Mach numbers approaching unity. 
Such Mach numbers characterise warm media, while colder 
parts of the ISM tend t o have supersonic velocities (see 
lElmegreen & Scald l2004h leading to even larger de nsity fluc- 
tuations jFalgarone et al.Lll998HPadoan et allfl997l) . . 

To a certain extent these problems can be circumvented by 
combining the density shift with a truncation of the residual 
negative tail. When we shift the density distribution e.g. by 
po = lo>, by adding this constant value, and discard all points 
falling below zero only 8 % of the points from the original dis- 
tribution are set to a zero value so that the statistical properties 
of the overall structure are hardly changed. In this way we can 
obtain positive densities and acr p /(p) ratio of about one, avoid- 
ing all problems from a dependency on the resolution and on 
the exact random numbers. One has to keep in mind, however, 
that the truncation of the density structure can have a notice- 
able influence on the scaling properties. The pure addition of 
the constant density does not affect affect the spectrum because 
it is scale-independent. 

To test the possible error introduced by the truncation of 
the distribution at a given density level we have analysed the 
truncated fBm's and compared them to the original spectra. 
The result is shown for a spectral index p - 2.6 and differ- 
ent truncation levels in Fig. |4] The actual shift of the density 
by po does not influence these spectra because the A-variance 
is insensitive to any constant offset. We see the pure truncation 
effect. For truncation levels of 2cr p and above the spectra are 



practically not changed. For truncation levels between 0.5cr and 
2cr the shape of the spectra is retained but they are shifted to 
lower absolute values. This can be explained by the reduction 
of the total variance in the data cubes which is visible in the 
scale-dependent A-variance as well. The original distribution 
was normalised to a variance of unity in this example whereas 
the truncation leads to reduced variances of 0.96, 0.76, and 0.56 
for the 2<x, lcr, and 0.5<x truncation levels, respectively. These 
are exactly the numbers by which the A-variance spectra in Fig. 
|4]are shifted relative to the original spectrum. Only for a trunca- 
tion level at 1 /4<x the slope of the spectrum is changed, i.e. the 
scaling behaviour of the structure is modified. In this case the 
absolute shift of the A-variance spectrum also does no longer 
exactly match the corresponding reduction of the total variance 
of the density distribution relative to the original value. 

Examining the resulting projected maps shows that the re- 
lation between the 3-D scaling and the 2-D scaling given in 
Eq. l|9} is also preserved down to truncation levels of 0.5cr. 
Corresponding studies for different spectral indices show that 
the A-variance spectra are least sensitive to truncations at low 
spectral indices, between 2 and 2.5, where even truncation lev- 
els of 0.25cr do not change the scaling behaviour and the re- 
lation between total variance and A-variance. At spectral in- 
dices close to four, in contrast, the 0.5<x truncation plot shows 
already significant deviations, so that we conclude that a negli- 
gible statistical impact on the scaling behaviour is only guaran- 
teed at truncation levels around lcr and above. The shift-and- 
truncate method to create positive densities is thus not perfect 
in terms of retaining the original scaling properties of the struc- 
ture, but the introduced deviations are still small compared to 
those introduced by the non-linear transformations discussed 
in Sect. 12.21 They would be hardly detectable in observed data, 
although we have to take them into account when performing a 
detailed quantitative analysis. 

Fig- E] shows two actual examples for the influence of the 
density zero level definition on the measured centroid velocity 
spectra. The scaling behaviour of the centroids was computed 
in terms of the A-variance spectra for three different shift-and- 
truncate levels of fBm generated density structures. To judge 
how far they reflect the original density or velocity structure, 
we have also plotted the A-variance spectra of these projected 
quantities multiplied with the mean square of the complemen- 
tary quantity to guarantee units equivalent to the centroids. 

The upper plot shows the combination of a shallow den- 
sity spectrum with a steep velocity spectrum, matching a situa- 
tion which is typically observed in molecular clouds (see Sect. 
12. 2> . The absolute shift of the curves for weighted centroids 
is mainly determined by the different values of (p 2 ) produced 
by different average densities. However, this shift does not in- 
fluence the characteristic scaling behaviour within the struc- 
ture. Looking at the slopes of the centroid spectra, we find a 
confirmation of the general considerations on the role of the 
density zero level po given above. If the density structure is 
dominated by a large average, i.e. in the case of po = 3cr p , the 
centroid velocities are basically given by a projection of the 
velocity structure, so that they reproduce the original velocity 
scaling behaviour. For lower average densities, i.e. a lower rel- 
ative contribution of the pure velocity projection given by the 
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Fig. 5. Comparison of the A-variance spectra of the weighted 
centroid velocities with the spectra of the original density and 
velocity structure for different truncation levels. To plot equiv- 
alent quantities the projected den sity is mu ltiplied by ( v 2 ) an d 
the projected velocity by (p 2 ) (see Esquivel & Lazarian, 2005). 
The (p 2 ) factor was computed for the 3cr-cut density cube. The 
corresponding plots for lcr and 0.25cr would be shifted down 
by a factor 5.2 and 13.5, respectively. The upper plot was com- 
puted from an fBm density structure with ft = 2.6 and a velocity 
structure with (3-3.1, the lower plot used the opposite spectral 
indices. 

first term in Eq. i ll Oi l, the centroid scaling becomes shallower 
with an exponent which is close to that of the velocity structure 
at large scales and an exponent which is close to that of the den- 
sity structure at very small scales and the lowest values of po. 
This plot seems to confirm the transition from purely velocity- 
dominated centroids to den sity-dominated cen troids as origi- 
nally interpreted bv lLazarian & Esquivellll2003h . 

If we consider, however, the opposite situation of spectral 
indices in the lower plot, we only find that the centroid scal- 
ing becomes less and less representative for the actual velocity 
structure when reducing the average density po. Their scaling 
does not tend towards the scaling of the column density struc- 
ture but becomes shallower as well. This fact is confirmed in 
all simulations with other combinations of spectral indices. At 
low values of o~ p /po, the centroids match the projected velocity 



Fig. 6. Same as Fig.|3]but for normalised centroids. They are 
rescaled by the factor (p 2 ) to obtain comparable dimensions. 

structure, whereas their scaling becomes shallower for lower 
average densities irrespective of the actual spectral index of the 
density structure. The spectral index of the density structure 
determines, however, at which scales the deviations occur. For 
density structures with a shallow spectral index, dominated by 
many small-scale fluctuations, the main effect occurs at small 
scales. In contrast, we find the main deviations at large scales, 
when the density structure has a steep index, representing a rel- 
ative dominance of large scale fluctuations. When interpreting 
changes in the slope of the A-variance spectrum of observed 
centroid maps we can thus use the known information on the 
projected density scaling to judge whether they represent an 
actual deviation of the velocity structure from self-similarity or 
whether they might be produced just by the centroid composi- 
tion effects. In general, we find a bigger impact on the overall 
centroid spectrum as shallower both spectra become. However, 
is not clear that the a shallo w velo city field is physically moti- 
vated (Es quivel & Lazarianl l2005). 

In Fig. [6] we show the same effects for the ordinary, nor- 
malised centroids instead of the weighted centroids used in Fig. 
|5] We find the same general behaviour as for the weighted cen- 
troids but differences in details. In all cases with high average 
densities, i.e. for cr p /po < 1 the normalised centroids provide 
a slightly better reproduction of the original velocity structure 
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than the weighted centroids. At lower densities, they are some- 
what weaker modified for shallow density spectra and some- 
what stronger modified for steep density spectra. The modifica- 
tions correspond approximately to the same effect that a change 
of o-p/po by a factor 1 .5 would have for the weighted centroids. 
In general we can state, however, that either both centroid def- 
initions reveal the true velocity structure or none of them. The 
direct retrieval of the velocity scaling from the A-variance spec- 
tra of the centroids will only succeed when the average cloud 
density is significantly larger than the density dispersion. In 
these cases the normalised centroids are marginally better than 
the weighted centroids. 

These results explain the differences and agreements be- 
tween the previous studies on velocity centroi ds discussed so 
far in t he l iterature. The studies of pVIiville-Deschenes et alJ 
(2003a) and L azarian & Esquivel ( 2003) used a relatively large 
average density and they indicated a good match between 
centroid spectra and projected velocity spectra. The mecha- 
nism of producing posi tive densities from fBm's by ad ding 
a large constant used bv lMiville-Deschenes et alJ (|2(303a) and 
Esa uivel et al.1 J2003h gave results that correspond to our results 
for applying the shift-and-truncate technique with a large aver- 
age density, i.e. when we add the 3cr density offset. Both cen- 
troid definitions follow the actual velocity scaling over a large 
range of scales, deviating at most at the very ends of the spec- 
tra in this case. One has to emphasise that this matching is only 
produced by adding a large po value, so that the general conclu- 
sion that c entroids are a good measure for th e velocity structure 
drawn by iMiville-Deschenes et alJ d2003al) and applied to in- 
terpret observational data bv lMiville-Deschenes et alJ (|2£)03b) 
does not hold for the general case of interstellar gas with sub- 
stantial density fluctuations. 

In cases with lower average densities, all using a combi- 
nation of steep velocity spectr a with shallow density s pectra , 
lOssenkopf & Mac Low) d2002l) . lLazarian & Esquivell (|2003), 
and Brunt & Mac Low! (120041) found centroid spectra which 
were shallower than the velocity spectrum. The hydrody- 
namic and magneto-hydrodyn amic turbulence models studied 
Ossenkopf & Mac Low (2002) were characterised by steep ve- 
locity spectra with /3 V ~ 4.0, shallow density spectra with 
P P * 2.5 . . .2.7, and a high density contrast with <x p /p > 5 
thus corresponding closely to the conditions for the low density 
curve in the upper plot of Fig. [6] With a limited dynamic range 
for fitting the A-variance spectra it is obvious that the flattening 
of the centroid spectra relative to the original velocity spectrum 
seen in the figure can be misinterpreted as a constant reduction 
of the slope by one. 

Our results can a lso explain the findings of 
iBrunt & Mac Low! ( 12004 studying the characteristics of 
velocity centroids of HD and MHD turbulence simulations as 
a function of Mach number. With the kn own relation betwee n 
Mach number and density dispersion (Padoa n et all LL997), 
their finding of a growing discrepancy between the average 
spectral index of the velocity distribution and of the centroid 
map with growing Mach number can be explained by the 
impact of an increasing cr p /po ratio, which reduces the relative 
contribution of the projection term. This is most clearly seen 
in the models of decaying turbulence where, e.g., for an initial 



o-p/po ratio of 1.1 the centroid spectrum is shallower by 0.8 
than the velocity spectrum whereas it is for the final <T p /p 
ratio of 0.5 only shallower by 0.1. We have to emphasise, 
however, that this approach cann ot explain the differences 
in the spectral indices obtained bv IBrunt & Mac Low! (|2004) 
for MHD models observed perpendicular or parallel to the 
main magnetic field direction. In these cases, the isotropy 
assumption used in our analysis is clearly violated. 



3.2. The velocity zero level 

The composition of weighted centroids (Eq.|5} is a priori sym- 
metric with respect to density and velocity. In the decomposi- 
tion in Eq. dlOi . we have assumed, however, that the velocity 
scale is chosen in such way that vo = while po > 0. To bet- 
ter understand the centroid behaviour it is useful to perform an 
experiment using velocity fields with vq > 0. From the symme- 
try of the problem, we expect that we find a centroid behaviour 
matching the density scaling for large average velocities vo in 
the same way as we find centroids matching the velocity scal- 
ing for large average densities po. Indeed, we obtain the new 
term vt (j dz 6p(x) xjdz 6p(x + t)) x in Eq. (HI if vo ^ 0. It 
contains the spectrum of the projected density fluctuations. In 
contrast to the density treatment, we do not apply any trunca- 
tion to the velocity structure when shifting it to vo > 0. The 
experiment thus provides an additional test for the significance 
of the truncation. If the simple shift of the velocity structure be- 
haves equivalent to the shift-and-truncate of the density struc- 
ture we can be sure that all effects result from the selection of 
the average values and not from the truncation. 

Fig.0shows the impact of different velocity offsets on the 
centroid A-variance spectra. A high average density, po = 3cr p 
was chosen, to guarantee that the centroid spectrum for vo = 
is dominated by the velocity structure. The curves for vo = 
are identical to the po = 3<x p curves in Fig. [5] When increas- 
ing the average velocity we find in the upper plot a transition 
to shallower spectra similar to the effect of a reduced average 
density in Fig. |5] The slope of the centroid spectrum remains 
close to the slope of the velocity spectrum at large scales and at 
small scales it takes the slope of the column density spectrum. 
In the lower panel we find as well that the centroid spectrum is 
more and more similar to the column density spectrum when 
increasing the average velocity. This is opposite to the effect of 
reducing the average density in the lower panel of Fig. [5] The 
adjustment of the average velocity reproduces the transition 
from velocity-dominated spe ctra to density- dominated spectra, 
as predicted by Lazarian & Esquivell (I2003I) . 

The equivalence of the impact of the velocity shift on the 
centroid spectra to the impact of the shift-and-truncate method 
for the density structures proves that the main change of the 
centroid spectrum is due to the added offsets and not due to 
the truncation of the density structure at its low density wing. 
Unfortunately, the numerical experiment cannot be exploited 
to derive the true velocity scaling when the average density 
is so small that the centroid spectra for vo = are "density- 
contaminated". By increasing vo we will only increase the con- 
tribution from the density scaling, which is already known from 
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Fig. 7. A-variance spectra of the centroid velocities computed 
using the density and velocity structures from Fig. [5] and dif- 
ferent v'o levels. A density structure with (p) = 3<x p was used 
here. 

the projected intensity maps, but we cannot remove the effect 
of the combination of density and velocity fluctuations. 

3.3. Further decomposition 

The results obtained so far show that the density zero level ba- 
sically changes the contribution from the first term in Eq. dlOl 
representing the pure projection of the velocity structure. In a 
next step we investigate the relative contribution of the other 
three terms to the deviation measured between the projected 
velocity spectra and the centroid spectra. The third and fourth 
term vanish if there is no cross-correlation between the density 
and velocity fields. For our independently generated fBm struc- 
tures this should be the case. We expect, however, that in every 
realization some accidental correlations occur so that the two 
terms are only negligible in the ensemble average. 

If the density field is known, we can obtain the second term, 
i.e. the combination of density and velocity fluctuations, by 
constructing an auxiliary density field p aux = p - po and com- 
puting the weighted centroids for this auxiliary quantity. As the 
average density of the auxiliary field vanishes, the derived cen- 
troids directly match the second term in Eq. dlOi . This proce- 



Fig. 8. Decomposition of the A-variance spectra of the 
weighted centroid velocities for the density and velocity struc- 
tures from Fig. fusing shift-and-truncate by 1/4 of the original 
standard deviation for the density structure. The solid line rep- 
resents the contribution from the velocity structure projected 
with a po weighting and the dash-dot line represents the cen- 
troid contribution from the auxiliary density fluctuation field. 
The sum of both terms (dash-dot-dot-dot) is very close to the 
measured centroid spectrum (dashed). 

dure is illustrated in Fig.[8]where we plot the centroid spectrum 
for the auxiliary field p aux , and compare the full centroid spec- 
trum obtained from the original density structure with the sum 
of this second term and the pure velocity scaling term. For the 
sake of comparison we also plot the spectrum of the projected 
density and velocity structure, where the velocity spectrum is 
multiplied here by pfj to represent exactly the first term in Eq. 
l[[0). The same combination of spectral indices as used in Fig. 
|5]was taken. A 0.25<x p shift-and-truncate level was used for the 
density structure, so that the centroid spectrum deviates con- 
siderably from the spectrum of the velocity fluctuations. 

For all studied combinations of spectral indices, we ob- 
tain a good match between the sum of the projected velocity 
spectrum and the centroid spectrum from the auxiliary field 
of density fluctuations with the full spectrum of the velocity 
centroids. Nevertheless, we find always a non-negligible dif- 
ference between the two curves, resulting from the acciden- 



12 



Ossenkopf et al.: The statistics of centroid velocities 



tal cross-correlations contributing to the third and fourth term, 
which are not contained in the sum. We also find that the second 
term, giving the combination of all fluctuations, has a spectrum 
which is always shallower than either of the projected spectra 
involved. For a steep spectrum of density fluctuations its slope 
turns even negative at large scales. This explains why the total 
spectrum of the weighted centroids is always shallower than the 
projected velocity spectrum, independent of the spectral index 
of the density spectrum. 

The computations have confirmed the theoretical expecta- 
tion, that the spectrum of velocity centroids consists of only 
two main contributions: the pure projection of the velocity 
structure determined by the average density po an d a shallow 
term mainly determined by the density fluctuations. In Sect. |4] 
we show how this decomposition can be exploited to measure 
the actual velocity structure from observed centroid maps if the 
cr p /po ratio can be estimated independently. 

3.4. Matching criteria 

We have seen that the single quantity giving the ratio between 
the strength of the density fluctuations and the average density 
Op/po is arj le to discriminate between the different behaviours 
of the centroid spectra. For low values of this ratio, the spec- 
tra are dominated by the actual velocity structure so that the 
3-D velocity scaling is preserved in observed centroid maps. 
For (Tp/po SS 0.5 the A-variance spectrum of the centroid map 
directly measures the spectral index of the underlying veloc- 
ity structure. For higher values, the centroid spectra are always 
shallower than the spectra from the projected velocity maps. 
They are produced by a combination of density and velocity 
fluctuations. 

In contrast to the suggestion of a density-dominated regime 
by ILazarian & Esauivell J2003I) . the systematic study of a 
wide range of combinations of spectral indices with the A- 
variance spectra shows no indications of a transition from ve- 
locity dominated centroids to density-dominated centroids, but 
rather a transition to "density-contaminated" spectra. Using 

a decomposition of structure functions similar to Eq. dlOt 

■ 1 1 1 I a 

Lazarian & Esquivel ( 2003) identified a term that indeed traces 
density fluctuations. They showed that in general centroids do 
not trace directly the velocity fluctuations. However, in their 
numerical tests, they use a combination of steep velocity and 
shallow density spectra, and disregard a cross term that is 
equivalent to the convolution of velocity and density fluctua- 
tions presented here. The shallow centroid spectrum was in- 
terpreted as the density spectrum. We have demonstrated that 
for the A-variance spectra a density-dominated regime arises 
only if we chose a velocity scale with an offset so that vq + 0. 
However, the combination of the facts that the centroid spec- 
trum is always shallower than the velocity spectrum and that 
most observed density spectra are shallower as well can give 
the false impression that centroids trace the density scaling for 
large ratios <r p po. 

Esqu ivel & Lazarianl J2005I) presented another criterion 
for a match between centroid and velocity scalings, namely, 
X 2 cr 2 , c » (v 2 )cr 2 ^ . They stated, however, that it is not clear how 



large the ratio X 2 cr 2 / ((v 2 )cr 2 ) should eventually be to guaran- 
tee that the centroids reliably represent the velocity structure. 
When applied to the overall data cubes we find that the distinc- 
tive power of the criterion is limited. In the examples plotted 
above we obtain for instance a ratio of 35 when using the 3<x 
shift-and-truncate level of the density distribution and a ratio 
of 3.0 for the 0.25<x shift-and-truncate level in the case of the 
shallow density and steep velocity spectrum. In contrast, we 
obtain corresponding ratios of 2.6 and 0.32, respectively, for 
the combination of steep density and shallow velocity spec- 
trum. In both cases the 3<x shift-and-truncate level gives a good 
match between centroid and velocity scaling while the 0.25cr 
level results in a very poor agreement. Thus the global crite- 
rion is poorly quantified. 

The criterion can be re written in a scale-dependent form 
(ILazarian & Esquivel 120031) : X 2 D Vc (l)/((v 2 }D h Jl)) » 1, 
when we consider the structure function of the two maps at 
a given lag I. One might assume that this criterion should hold 
as well for A-variance spectra because of their similar scaling 
properties. Then a ratio X 2 a\ v (l)j((y 2 )(j\ l (/)) much larger 
than unity indicates a good match of the centroid A-variance 
spectrum with the true velocity spectrum. The denominator 
grows relative to the numerator with increasing scales when 
the density spectrum is steeper than the velocity spectrum. In 
this case, matched by the lower panels of Figs. 15181 the largest 
deviations of the centroid spectrum from the velocity spec- 
trum should occur at large scales, whereas the slopes of the 
A-variance spectra should match at small scales. This is indeed 
the behaviour that we observe in these figures. For the oppo- 
site relation of spectral indices, where the velocity spectrum is 
steeper than the density spectrum, as seen in the upper panels 
of the figures, the ratio is growing towards larger scales, and in 
fact we find the best matches of the scaling behaviour at large 
scales and the main deviation at small scales. 

Using the A-variance spectra in Fig. EI we can evaluate the 
criterion by eye from the plots. When the curves for the cen- 
troids fall well above the dotted line giving the density spec- 
trum, the centroid spectrum should be a reliable tracer of the 
velocity structure. The same test can be performed in the anal- 
ysis of observed data, because the A-variance spectra of the in- 
tensity and the centroid velocity maps and the average velocity 
dispersion are easily measured in observed line data. However, 
we find that the actual significance is also limited. In the up- 
per panel of Fig.[5]with the shallow density and steep velocity 
spectrum, we find that the velocity spectrum is reproduced by 
the centroid spectrum when the centroid A-variance exceeds 
the values from the column density structure by about a factor 
four, whereas for the steep density and shallow velocity struc- 
ture we get a good match even if the centroid curve falls just 
above the column density spectrum. For other combinations of 
spectral indices we find that a ratio of two is sufficient to guar- 
antee a match between centroid spectrum and velocity spec- 
trum as long as the density spectrum is very steep (J3 p > 3.5), 
whereas ratios as high as 100 may be required to guarantee a 
match when the density spectrum has an index shallower than 
2.5. When using the normalised centroids in Fig.|6]we cannot 
derive an equivalent criterion to estimate the match between 
centroid scaling and velocity scaling based on the measured 
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map spectra. This is a clear practical advantage of the weighted 
centroids. 

Thus we can basically confirm the criterion, when ap- 
plied in its scale dependent form to A-variance spectra, but 
have to emphasise that there is no single value for the ratio 
where the transition between velocity-dominated and "density- 
contaminated" behaviour appears, but that the exact shape of 
the density spectrum has to be taken into account. 

3.5. Comparison of A-variance spectra and 
structure function 

As the structure function is related to the autocorrelation func- 
tion, the decomposition in Eq. (II Oi applies as well to the contri - 
butions to the structure function. Laz arian & Esauivell (12003 ). 
however, have shown that the second term representing the 
combination of density and velocity fluctuations can be further 
split into two separate contributions in terms of the structure 
function. As one of them represents the pure density fluctua- 
tions, they suggested that the structure function can undergo 
a transition from a velocity-dominated spectrum to a density- 
dominated spectrum. 

To test this behaviour we have repeated the experiments 
shown in Figs. 13 to for structure functions. In general we 
expect to see clear deviations from power-laws as the projec- 
tion of structure functions results always in broken power laws. 
Structure functions of 2-D projections can be represented by 
two asymptotic power laws: one at small lags (|/| «: Ztot) hav- 
ing a spectral index f3 - 2 for both shallow and steep spectra, 
and another one at large lags (|/| » Ztot) with a spectral in- 
dex H — 3 for steep spectra an d (constant) for shallow spectra 
(Esa uivel & Lazarian 120051 . Taking the general limitation of 
a restricted dynamic range of scales, both in the fBm simula- 
tions and in most observed maps, the 2-D structure functions 
will always fall in the transition between the two asymptotes 
so that their slope cannot be reconciled directly, preventing a 
direct recovery of the underlying 3-D statistics. Hence, no sim- 
ple inversion of the projection problem is possible. Compared 
to the A-variance spectra, the spectra of structure functions are 
thus always somewhat more curved with steeper slopes at small 
lags and shallower slopes at large lags, but in spite of the dif- 
ferent analytic decomposition of the structure function ofcen- 
troid velocities demonstrated by Esauivel & Lazarian (2005), 
the general behaviour is always very similar to the A-variance 
spectra. 

The measured changes with respect to variations of the den- 
sity and velocity zero level is also almost identical to the be- 
haviour shown in Figs.|5]toQ For the centroids obtained from 
the density structure with the 3cr p shift-and-truncate level we 
find a very good match of the structure functions of the cen- 
troids and the projected velocity structure. If the density disper- 
sion, however, is in the order of the average density the spectra 
flatten with main deviations at large scales when the density 
spectrum is steep and at small scales when it is shallow. When 
comparing A-variance spectra and structure functions in detail, 
we find that the centroid structure functions resemble the true 
velocity structure always slightly better than the A-variance 



spectra. This might be partially due to the somewhat lower sen- 
sitivity of the structure function to cha nges in the power spec- 
trum a t particular scales as found by lOssenkopf & Mac Low! 
(2002k but might indicate also a slight advantage of the struc- 
ture function compared to the A-variance spectra when applied 
to centroid maps. 

There is again no transition from a velocity-matching be- 
haviour to a density-matching behaviour, but rather a density- 
contaminated structure with a spectrum which is shallower 
than the true velocity spectrum. We fin d as well a confirma- 
tion o f the scale dependent criterion of Esauivel & Laza rian 
d2005l) for a match between centroid spectrum and true veloc- 
ity spectrum. Here, the critical ratio X 2 D Vc (l)/((v 2 )Di iM (l)) for 
a match between centroid structure function and velocity struc- 
ture function for a particular combination of spectral indices is 
always somewhat smaller compared to the A-variance spectra. 
For steep density spectra, a ratio of one seems to be always suf- 
ficient, whereas for shallow density spectra a ratio of 20 may be 
required to guarantee a velocity-dominated centroid behaviour. 

4. The derivation of the velocity structure from 
density-contaminated centroids 

Whenever the average density of the medium is too small so 
that the centroid spectrum does no longer reflect the underlying 
velocity spectrum, we can deduce the true velocity spectrum 
from measured centroids only when we find a way to compute 
the second term in Eq. JlQi and when the last two terms pro- 
duced by the accidental correlations are negligible. Based on 
the results of the decomposition shown in Fig. [8] we propose 
an iteration scheme which computes the second term, i.e. the 
convolution of the two fluctuation spectra, from the first term, 
i.e. the velocity projection weighted with pi, obtained in a pre- 
vious iteration step, neglecting the small contribution from the 
other two terms. 

As a first step to compute the fluctuation term a three- 
dimensional fluctuation structure has to be constructed which 
matches the scaling behaviour of the measured column den- 
sity structure. This can be done in the following way: From the 
measured A-variance spectrum of the density projection, i.e. 
the spectrum of the intensity map, we can compute the three- 
dimensional A-variance spectrum by de-projecting it according 
to the results from Sect. 12.51 By translating this spectrum by 
k - Zm/l mt0 a spherically symmetric power spectrum P(k), 
we create a new fBm-like structure using this power spectrum 
and random phases. This new structure should match the scal- 
ing behaviour of the input 2-D A-variance spectrum. 

This is illustrated in Fig.|9] We demonstrate the quality of 
this construction for two examples. In the first case we use the 
A-variance spectrum of the projection of a known fBm with 
a spectral index /? = 2.6, in the second case we start from 
an artificial spectrum given by a power law corresponding to 
/3 = 3.7 up to 17 pixels and an V 1 decay above. In both cases 
we create the corresponding 3-D fluctuation field, compute its 
2-D projection and the A-variance spectrum of the projection. 
Comparing this derived spectrum with the input spectrum in 
Fig-IUgives an impression of the quality of the reconstruction. 
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Fig. 9. Two examples for the construction of fluctuation fields 
from a projected A-variance spectrum. The figure compares the 
input spectra with the spectra obtained from the A-variance 
analysis of the projection of the computed fluctuation fields. 

For the power-law input spectrum we find an almost per- 
fect match with small deviations due to the artificial gridding 
of the fluctuation field, numerical uncertainties, and statistical 
fluctuations. In contrast, the spectrum composed of two power- 
laws is less accurately reproduced. The A-variance spectrum 
of the fluctuation field shows a broader peak and approaches 
the original spectrum only at lags relatively far apart from the 
peak. This broadening is due to the convolution of the power 
spectrum with the filter function in Eq. which was ignored 
in the simple translation of the A-variance spectrum back into 
a power spectrum described above. In principle we could try to 
include a corresponding deconvolution to make the approach 
fully self-consistent, but the reasonable agreement between the 
two curves even in this extreme case shows that this additional 
refinement is not needed. The example was chosen to be ex- 
treme in the sense, that we have a sharp turn from a steeply 
increasing spectrum into the steep decay of the A-variance rep- 
resenting completely uncorrelated structures. In all cases with 
wider peaks, the agreement between the original spectrum and 
the derived spectrum is better, although the general tendency 
remains that the peak in the derived fluctuation spectrum is al- 
ways slightly too broad. The actual quality of the construction 
of the fluctuation field from the A-variance spectrum will thus 
fall between the two extremes shown in Fig. [9] 

The fluctuation field constructed in this way has a zero av- 
erage, so that we can use it directly as the auxiliary field to 
compute the centroids for p - po in Fig. [8] i.e. the second 
term in Eq. JlOb . Unfortunately, the unknown field of veloc- 
ity fluctuation enters as well into this term so that an iteration 
scheme is required: we start from the measured centroid spec- 
trum, assuming that it is purely determined by the projection 
of the velocity field, divide by pS? and construct a fluctuation 
field for the velocities in the same way as described above for 
the density fluctuation field. From the convolution of the two 
fluctuation fields we estimate the A-variance spectrum of the 
correction term. Subtracting this spectrum from the measured 
centroid spectrum then provides the next estimate for the pure 
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Fig. 10. A-variance spectra of the weighted centroid veloci- 
ties for the density and velocity structures from Fig. [8] (dotted 
lines). The solid line represents the projected velocity structure 
with a po weighting. The dashed and the dash-dot lines rep- 
resent the correction term and the derived projected velocity 
contribution at the end of the iteration. In the ideal case this 
converged solution should agree with the spectrum from the 
original velocity structure. 

projection of the velocity field. This can be used again to de- 
termine the pure fluctuation term in the centroids and so on. 
The iteration is stopped when the velocity spectrum obtained 
in subsequent steps remains constant within 1 % 7 

An example for the result of this iteration is displayed 
in Fig. [K)] for the centroid spectra obtained from the combi- 
nation of shallow density and steep velocity fields and vice 
versa as shown in Fig. [8] The general recovery of the pro- 
jected velocity structure is quite satisfactory. The absolute mag- 
nitude of the fluctuations is, however, somewhat too small, and 
for the combination of shallow density with a steep velocity 
spectrum the derived overall velocity spectrum is also slightly 
steeper than the original spectrum. These remaining deviations 
should stem from the accidental correlations between density 



7 The exact value of the convergence criterion is not important, it 
only changes the number of required iterations. We found that the re- 
sults obtained for smaller error limits cannot be distinguished by eye 
from the 1 % limit results. 
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Fig. 11. Derivation of the projected velocity spectrum from a 
measured centroid spectrum when applying a 20 % variation to 
the density po used in the iteration scheme. 

and velocity field, expressed in the higher terms of Eq. ( I10i . 
Altogether, the iteration scheme has proven to be a reliable 
method to recover the original velocity spectrum from a mea- 
sured centroid spectrum, when the projected density structure 
and the average density are known. In all fBm combinations 
tested here, the overall slope of the derived velocity spectrum 
agrees with that of the original spectrum within 0.1. This is 
sufficient to distinguish between different turbulence models 
felmegreen & ScalolEooH) . 

A major drawback of the method is the need for an accurate 
estimate of the average density in the considered interstellar 
cloud. This is not easy to obtain from the projected density in 
observational data because the line-of-sight extent of a cloud is 
often not known. This can be overcome in clouds with a known 
geometry or by excitation studies of molecular tracers sensitive 
to particular densities. However, very accurate estimates will 
be always difficult. 

Thus we have studied the influence of an error in the de- 
termination of the average density on the reconstruction of the 
velocity structure. In Fig.^2 we nave repeated the experiment 
shown in the upper plot of Fig. ^3 when increasing and de- 
creasing the used average density relative to the actual value by 
20 %. The result shows the same tendencies discussed in Sect. 
O When the average density is overestimated, the centroids 
are thought to better resemble the scaling of the underlying ve- 
locity structure. The velocity fluctuations are underestimated 
because they are obtained by dividing the spectrum by a p 2 , 
value which is too large. The computed correction term is too 
small and the derived velocity spectrum falls above the actual 
spectrum and is too shallow. If the average density is underesti- 
mated, we correct the centroid spectrum with an overestimated 
fluctuation spectrum, so that the derived velocity spectrum is 
too small and too steep. For steep velocity spectra a change of 
the average density by 20 % corresponds to a change of the av- 
erage exponent of the spectrum by 0.25. For shallow spectra, 
the influence is somewhat smaller. 

Thus we can conclude that it is possible to the retrieve the 
actual velocity spectrum from measured centroids even if the 



centroid spectrum is density-contaminated, but the accuracy of 
this retrieval depends critically on an knowledge of the average 
density in the cloud. Independent measures of the <r p /po ra- 
tio are required. Methods to accurately derive the density from 
multi-line observations ha ve been s u ccessfully dev elo ped an d 
applied bvlSchrever et al](e.g.ll997l) : lRichter et ai1 re.g. l2003h : 
ISonnentrucker et alJ (e.g. 120031) . Thev are based on the com- 
bination of information from different species tracing a wide 
range of critical densities, but the accurate determination of po 
still remains a challenging task. 

5. Conclusions 

We have shown that the A-variance analysis is an appropriate 
tool to characterise the scaling properties of both velocity cen- 
troid maps and the underlying three-dimensional velocity field. 
By directly reflecting the power spectrum of fluctuations and 
preserving a power-law behaviour through the projection the 
A-variance is well suited to quantify the properties of interstel- 
lar velocity fields. The fact that velocity centroids may not re- 
flect the velocity statistics was always a concern for turbulence 
research. The disagreement between the aforementioned statis- 
tics was_dj^cu^se^_ajready_by [Ossenkgjj^^Mi^Ls^ |2002) 
and lBrunt & Mac Low! J2004I) '). We successfully tested the cri- 
terion for the validity of centroids as measures of velocity 
statistics suggested by Lazar ian'& Esquiveilll2003l) . 

We find that the most accurate criterion determining 
whether a centroid spectrum reflects the velocity scaling prop- 
erties is a small ratio between the density dispersion and the 
mean density. For values below 0.5 the centroid spectra match 
the underlying velocity structure. Here, the centroids are deter- 
mined by the pure projection of the velocity field. At higher 
cr p /po ratios the mutual convolution of density and velocity 
fluctuation contributes a main term. Based on this knowledge, 
we can qualitatively explain all the differences in the interpre- 
tation of centroid spectra found in the literature. 

Without knowing the average density in the considered 
medium we can test whether a centroid spectrum reflects true 
veloci ty structure using the criterion bv lLazarian & Esauivell 
( 2003J) that XV 2 (Z) » <v 2 )cr 2 (/) when the centroid spectrum 
is velocity-dominated. Although derived for the structure func- 
tion it holds for the A-variance as well. However, there is no 
single value by how much the left hand side has to exceed the 
right hand side. We have confirmed the criterion by numeri- 
cal experiments and found that factors above two are sufficient 
in case of steep density spectra but factors up to 100 may be 
required for extremely shallow density spectra. Currently, ob- 
servations and simulations of interstellar turbulence show that 
both steep and shallow regimes may occur with density spectral 
indices ranging from about 2.5 to 3.3 (see Sect. 12. 2> . 

We do not see a transition from velocity-dominated to 
density-dominated spectra at lower densities, but rather a tran- 
sition to "density-contaminated" spectra which are systemati- 
cally shallower. The flattening of the centroid spectra relative 
to the true velocity structure in the general case can be easily 
misinterpreted as a transition from a velocity-resembling to a 
density-resembling spectrum because in interstellar turbulence 
the density spectra are often shallower than the velocity spec- 
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tra. In any case, the density structure can be obtained directly 
from column density maps. A density-dominated spectrum oc- 
curs only if the velocity scale was chosen in an unfortunate way 
so that the average velocity is not negligible relative to the ve- 
locity dispersion. By adjusting the velocity frame such that the 
average line is centred at zero, this term can always be elimi- 
nated. 

Whenever the centroid spectrum is velocity-dominated, the 
A-variance analysis is a simple and robust tool to directly in- 
fer the velocity scaling from the centroid map. The exponent 
of the A-variance spectrum is the exponent of the power spec- 
trum of the velocity fluctuations reduced by two. Although, the 
second order structure function is connected to the power spec- 
trum by a different functional behaviour, we find a very similar 
behaviour when applied to centroid velocities. All general con- 
clusions apply there as well. However, the structure functions 
of projections of power-law power spectra are always curved, 
so that a direct fit of the exponent is more difficult. Moreover, 
we find that, although our analytical decomposition of the ve- 
locity centroids is only valid for weighted centroids, the nor- 
malised centroids behave qualitatively in the same way so that 
they can be used as well to derive the velocity structure when- 
ever the centroid spectrum is velocity-dominated. 

We provide an iteration scheme to derive the actual velocity 
structure from the centroid maps in all cases where the aver- 
age density is known, even if the map is density-contaminated. 
An accurate determination of the power spectrum of the ve- 
locity fluctuations depends on three conditions: the correlation 
between density and velocity structure can be neglected, the 
dynamic range of length scales covered by the map is sufficient 
to compensate for statistical fluctuations at particular lags, and 
the average density po can be estimated with a high accuracy. 

A different iteration scheme can be developed using the 
structure function instead of the A-variance. In this cas e the de- 
composition proposed by lEsquivel & Lazarianl J2005t) can be 
used to obtain a scheme which is less sensitive to the knowl- 
edge of the average density, but a considerably more complex 
approach is needed to evaluate the projection effects. This will 
be the topic of a subsequent paper. 

Acknowledgements. VO was supported by the Deutsche Forschungs- 
gemeinschaft through grant 494A. AE acknowledges support from 
the NSF grant AST-0307869 and the Center for Magnetic Self- 
Organization in Laboratory and Astrophysical Plasmas, and Mexico's 
Consejo Nacional de Ciencia y Tecnologfa. AL is supported by NSF 
grant AST0307869. We thank an anonomous referee for comments 
helping to lay out many aspects in a clearer and more precise way. 
We have made use of NASA's Astrophysics Data System Abstract 
Service. 

References 

Bensch R, Stutzki J., & Ossenkopf V. 2001, A&A 366, 636 
Brunt CM., Heyer M.H. 2002, ApJ 566, 289 
Brunt CM., Mac Low M.-M. 2004, ApJ 604, 196 
Chappell D., Scalo J. 1999, MNRAS 310, 1 
Cho J., Lazarian A. 2003, MNRAS 345, 325 
Cho J. & Lazarian A. 2005, Theor. and Comput. Fluid 
Dynamics 19, 127 



Combes F. 2000, in: The Chaotic Universe, ed. V.G. 

Gurzadyan, R. Ruffini, World Sci, 143 
Deshpande A.A., Dwarakanath K.S., Goss W.M. 2000, ApJ 

543, 227 

Elemgreen B.G, Kim S. Staveley-Smifh L. 2001, ApJ 548, 749 
Elmegreen B.G., Scalo J. 2004, ARAA 42, 21 1 
Esquivel A., Lazarian A., 2005, ApJ , 631, 320 
Esquivel A., Lazarian A., Pogosyan D., Cho J. 2003, MNRAS 
342, 325 

Falgarone E., Pineau Des Forets G, & Roueff E. 1995, A&A 
300, 870 

Falgarone E., Panis J.-E, Heithausen A. et al. 1998, A&A 331, 
669 

Falgarone E., Hily-Blant P., Levrier F. 2004, Ap&SS 292, 89 
Goldreich P., Shridhar S. 1995, ApJ , 438, 763 
Higdon J.C 1984, ApJ 285, 109 

Huber D. 2002, Nonequilibrium, Self-Gravity and Fragmented 

Interstellar Medium, PhD thesis No. 3348, Univ. Geneve 
Jenkins E. 2004, Ap&SS 289, 215 
Kleiner S.C. & Dickman R. L. 1985, ApJ 295, 466 
Kim J., Ryu, D. 2005, ApJ 630, L45 

Langer W.D., Wilson RW. & Anderson C.H., 1993, ApJ 408, 
L45 

Lazarian A. 1995, A&A , 293, 507 
Lazarian A., 2004, JKAS 37, 563 
Lazarian A., Esquivel A. 2003, ApJ 592, L37 
Lazarian A., Pogosyan D. 2000, ApJ , 537, 720 
Levrier, F. 2004, A&A 42 1 , 387 

Lis D. C, Pety J., Phillips T. G, & Falgarone E. 1996, ApJ 463, 
623 

Mac Low M.-M. & Ossenkopf V. 2000, A&A 353, 339 
Miesch M.S., Bally J. 1994, ApJ 429, 645 
Miesch M. S. & Scalo J. M. 1995, ApJ 450, L27 
Miesch M. S., Scalo J., & Bally J. 1999, ApJ 524, 895 
Miville-Deschenes M.-A., Levrier F, Falgarone E. 2003, ApJ 
593,831 

Miville-Deschenes M.-A., Joncas G, Falgarone E., Boulanger 

F, A&A 411, 109 

Munch G. 1958, Rev. Mod. Phys. 30, 1035 

Ossenkopf V. 2002, A&A 391, 295 

Ossenkopf V. & Mac Low M.-M. 2002, A&A 390, 307 

Ossenkopf V., Bensch F, & Stutzki J. 2000, in: Gurzadyan V. 

G. & Ruffini R. (eds.), The Chaotic Universe, World Sci., 
p.394 

Ossenkopf V., Krips M., Stutzki J., 2005, A&A , submitted 
Padoan P., Jones J.T., & Nordlund A.P 1997, ApJ 474, 730 
Padoan P., Boldyrev S., Langer W. & Nordlund A.P. 2003, ApJ 
583, 308 

Peitgen H.-O. & Saupe D., 77ie Science of Fractal Images, 

Springer, Berlin, 1988 
Richter P., Sembach K.R., Howk J.C. 2003, A&A 405, 1013 
Schreyer K, Helmich FP, van Dishoeck E.F., Henning T. 

1997, A&A 326, 347 
SonnentruckerP, Friedman S.D., Welty D.E., YorkD.G, Snow 

T.P, 2003, ApJ 596, 350 
Stanimirovic S., Lazarian A. 2001, ApJ 551, L53 
Stenholm L.G., 1984, A&A 137, 133 



Ossenkopf et al.: The statistics of centroid velocities 



17 



Stutzki J., Bensch F., Heithausen A., Ossenkopf V., & Zielinsky 

M. 1998, A&A 336, 697 
Zank G.P., Matthaeus W.H. 1992, J. Plasma Phys. 48, 85 



